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Abstract 

In order to study the origin of high-frequency quasi-periodic oscillations observed 
in X-ray binaries, Kato (2004) suggested a resonant excitation mechanism of disk 
oscillations in deformed disks. In this paper, we investigate numerically, following his 
formulation, whether trapped g-mode oscillations in a warped disk, where the warp 
amplitude varies with radius, can be excited by this mechanism. For simplicity, we 
adopt Newtonian hydrodynamic equations with relativistic expressions for the char- 
acteristic frequencies of disks. We also assume that the accretion disk is isothermal. 
We find that the fundamental modes of trapped g-mode oscillations with eigenfre- 
quencies close to the maximum of epycyclic frequency are excited. The intermediate 
oscillations found are isolated in a narrow region around the resonance radius. After 
varying some parameters, we find that the growth rate increases as the warp ampli- 
tude or the black hole spin parameter increases, while it decreases as the sound speed 
increases. 
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1. Introduction 

The Rossi X-ray Timing Explorer (RXTE) satellite has detected high-frequency quasi- 
periodic oscillations (QPOs) in X-ray binaries. They are kilohertz (kHz) and hectohertz (hHz) 
QPOs in neutron-star low-mass X-ray binaries, and high-frequency (HF) QPOs (> 100 Hz) in 



1 



black-hole X-ray binaries. HF QPOs occur at fixed frequencies and the appearance is correlated 
to the state of the sources, that is, they appear only in high-luminosity states where L > O.ILe, 
with Le being the Eddington luminosity. They are suggested to be phenomena originating in 
a strong gravitational field, which occur in the innermost region of the relativistic disk. HF 
QPOs are regarded as being a powerful tool to explore the mass and spin of the central black 
hole, and also to explore the physical states of the innermost region of the accretion disk. 

Various models have been proposed to explain HF QPOs. For example, Abramowicz & 
Kluzniak (2001) and Kluzniak & Abramowicz (2001) pointed out the importance of resonant 
processes as being the cause of HF QPOs. In their model, HF QPOs are the result of resonant 
couplings between the vertical and horizontal epicyclic oscillations at a particular radius. There 
is, however, uncertainty concerning the excitation process of the oscillation system as a whole. 

Kato (2004; 2008a; 2008b) proposed a model where HF QPOs are regarded to be disk 
oscillations resonantly excited in deformed disks. The deformation considered is a warp or an 
eccentric disk deformation in the equatorial plane. An outline of the model is as follows. A 
non-linear coupling between a disk oscillation (hereafter, an original oscillation) and a deformed 
part of the disk (warp or eccentric deformation) causes some forced disk oscillations (hereafter, 
intermediate oscillations). The intermediate oscillations make a resonant coupling with the disk, 
and then feedback to the original oscillation. Since the nonlinear feedback process involves a 
resonance, the original oscillation is amplified or damped. There are two kinds of resonances, 
i.e., horizontal resonance (Lindblad resonance) and vertical resonance. When the deformation 
is a warp, the resonance that can excite p- and g-mode oscillations is the horizontal resonance 
(Kato 2004; 2008a; 2008b). Hence, in this paper we consider the horizontal resonance alone. In 
the Keplerian disks this resonant process works only when the disk is general relativistic. That 
is, a non-monotonic radial distribution of epicyclic frequency is necessary for the appearance 
of resonance and for trapping of oscillations. Recently, Ferreira & Ogilvie (2008) generalized 
Kato's idea on the excitation of trapped g-mode oscillations and made detailed numerical 
calculations of the growth rates of the oscillations in the case of rotating black holes. Their 
results indicate that the coupling mechanism can provide an efficient excitation of trapped g- 
mode oscillations, provided that the global deformations reach the inner part of the disk with 
non-negligible amplitude. 

In this paper, we consider nearly the same problem examined by Ferreira & Ogilvie 
(2008), following the formulation by Kato (2008a; 2008b). The purpose of this paper is to 
make comparison of results obtained from two different approaches, Ferreira & Ogilvie's fully 
numerical approach and Kato's analytical one. We especially consider the case where the disk is 
warped. In section 2, we summarize the basic equations for original oscillations and intermediate 
oscillations, and an analytical expression for growth rate, which we use in later sections. In 
section 3, we present our numerical results, and section 4 is devoted to conclusions and a 
discussion. A detailed formulation of the nonlinear coupling terms is provided in Appendix. 
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2. Summary of Basic Equations 



2.1. Original Oscillations 

We assume an unperturbed disk to be steady and axisymmetric with a flow Uq. Then, 
by using a displacement vector a weakly nonlinear hydrodynamical equation for adiabatic 
perturbations is written as (Lynden-Bell and Ostriker 1967) 

p ^| + 2p (u • V) ^ + L(0 = Po C(Z, 0, (1) 
where L(£) is a linear Hermitian operator with respect to given by 

= Po («o • V) (u • V) £ + po (£ • V) (V^ ) + V [(1 - ri)podiv£] 

-p V(div^) - V [(£ • V)p„] + (£ • V) (Vpo), (2) 

where Po( r ) an d Po( r ) are the density and pressure in the unperturbed state, respectively, T 1 is 
the barotropic index specifying the linear part of the relation between Lagrangian variations 5p 
and 5p, and ipo is a general-relativistic gravitational potential. The right-hand side of equation 
(1) represents the weakly nonlinear terms. For F± = 1 (isothermal perturbations), it is given by 
(Kato 2004) 

€) = -^^V*,) - A ( Po i|v & ) . (3) 

We write the displacement vector £ in the form of a normal mode, 

€(r,t) = £(r,z)eiq)[i(ut-m<p)], (4) 

where u is the frequency and m is the azimuthal wavenumber. In the case where the disk 
is isothermal in the vertical direction and the perturbations have a short radial wavelength 
compared with the radial characteristic length of the unperturbed disks, we can separate £ into 
r-, ip-, and z-components as (Okazaki et al. 1987) 

Ur,z)=^ n (r)H n (z/H), (5) 

i v (r,z)=^ )n (r)H n (z/H), (6) 

Ur^)=L, n {r)H n ^{z/H), (7) 

where n is a non- negative integer characterizing the number of node(s) of the oscillation in 
the vertical direction and 1-L n is the Hermite Polynomial of argument z/H, with H being the 
vertical scale-height of the disk, which is given by H = c s /Q±, where c s is the sound speed and 
Q± is the vertical epicyclic frequency. 

We express the r-, cp-, and z- components of the homogeneous parts of equation (1) as 
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£ r>n -i2n(w-mn)£ v , n + tfj_H-^ = 0, (8) 



(cj — mfl) 2 ^^ + i2Q(u — mf])£ r ,ji = 0, (9) 
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Fig. 1. Warp function W(r) for c s = 10 3 , Wo = 0.01, and a = 0.5. 
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Note that these equations are the same as equations (23) - (25) of Kato (2008b). From these 
equations we have the basic equation describing the original oscillation (Okazaki et al. 1987): 

d 2 lr,n [(w — mVt) 2 — k 2 } \nVt\_ — [uj — mf2) 21 



(11) 



dr 2 Q 2 ± H 2 (co-mn) 2 
where Q is the (Keplerian) angular frequency of disk rotation, and k is the horizontal epicyclic 
frequency, respectively In the general relativistic potential around a rotating black hole, these 
frequencies are given by (e.g., Kato 1990) 



K = Q 
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(12) 
(13) 
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where a (—1 < a < 1) is the central star spin parameter, and r is in units of the gravitational 
radius R g , where R g = GM/c 2 with M being the mass of the central star. 

2.2. Warp and Intermediate Oscillations 

As mentioned in section 1, our resonant excitation model for the high-frequency QPOs 
assumes the disk to be warped or to have an eccentric deformation. In this paper, we con- 
sider the case where the disk is warped. The Lagrangian displacement associated with the 
deformation, £ w (r,t), is given by 



£ (r,t) =exp(-im w <p)£ 



w 



[r,z 



(15) 



where m w = 1. The warp solution is described by u\f = — W(r)Qz, uZ = iW(r)zd(rQ)/dr, and 
= W(r)Qr, where , ujf , and are the Eulerian perturbation of the velocity field and 
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W(r) is the solution of 
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(Papaloizou & Lin 1995; see also Ferreira & Ogilvie 2008). As the boundary condition, we 
specify the value of warp amplitude at the inner radius, Wo, and set dW/dr = there, following 
Ferreira & Ogilvie (2008). Here and hereafter, we take the inner radius to be the radius of 
the Innermost Stable Circular Orbit (ISCO), which is the radius where the horizontal epicyclic 
frequency, k, goes to zero. Figure 1 shows a typical warp function W(r), where c s = 10 -3 , 
W = 0.01, and a = 0.5. 

* w 

Using the warp function W(r), we can write £ as 

£w = -iW(r)z = -iW(r)HH!(z/H), (17) 

$ = -W(r)z = -W{r)HUi{z/H), (18) 

|f = iW(r)r = iW(r)rH (z/H). (19) 

Nonlinear coupling between the original oscillation, characterized by (u,m) and the 
deformation £ w characterized by (0,1) introduces intermediate oscillations, £ mt , which can be 
expressed as 



int , 



r,t) = exp [i(ut — mip)] £ 



int 



(20) 
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where m = m + 1 or m — 1. Here £ + represents the intermediate oscillations with m = m 

a W * int 

resulting from the coupling between £ and £ , while represents those with fh — m — 1 

- w* 

resulting from the coupling between £ and £ , where the asterix represents the complex 

~ int 

conjugate. To consider these various coupling cases separately, we write £ ± in the following 
form: 

tX(r,z)=tX,n{r)Kn{z/H), (21) 

(22) 

tX(r,z)=tXA r )Kr~Uz/H), (23) 

where n = n + 1 or n — 1 for the coupling with the warp. 

The nonlinear coupling terms are separated into terms proportional to exp[i(ojt — rhtp)] 
and T-Lfi{z/H). In the case of coupling through £ w , we write the coupling terms as 
1 
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C(£ w , £)\ r = po E A-,+,n(r) exp [i{ut - m<p)} H^z/H) + • • • , (24) 

h 

C(£ w ,£j\ =p Q J2 A v,+dr)eM^t-™<P)]nn(z/H) + ---, (25) 
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C(t £ w ) + C(£ w , £)] = Po E A Zt+ ,rSr) exp [i{ut - rh V )} H&-1 (z/H) + • • • , (26) 



where m = m + 1 and H denotes terms orthogonal to both Hfi and Jin-i- The subscript 

+ of A's denotes that they are related to the yj-dependence of exp[— i(m+ 1)^]. Note that 
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these equations are the same as equations (29) - (31) of Kato (2008b). In the case of coupling 
through £ w *, the nonlinear coupling terms, (l/2)p [C($ ) ,£ w *) + C(£ H/ *,£)], can be expressed 
in similar forms, by introducing Ar-^A^,-^, and ^4 z ,-,n instead. Then, the basic equations 
describing intermediate oscillations are written as 

d 2 



-(u-m{iy + K 2 -4n 



2 c 2 



'dr' 2 



*int 



£X ft -i2fi(o,-mfi)^ ±ifi 
+ tfH^^ = A 



dr 



dP nt 
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(27) 
(28) 

(29) 



where both cases of m = m + 1 and rh — m—1 are written together. 

We now eliminate ^±,h anc ^ C n ±,n fr° m the set of equations (27) - (29) to derive an 
equation with respect to £* n ± h . By using the approximation that the radial wavelength of 
oscillations, A, is shorter than r, i.e., A r, which has been used in deriving equations (8) - 
(10), we have 



(u-rhn) 2 + K 2 ]^l, 
2VL 



c 2 (u-mn) 2 d 2 (^%; 



(uj — mtty 



A 



r,±,ra 



-.4 



uj — mVL {p ' ±,n 1 (uj — mVL) 2 — n{l'j_ dr 
The horizontal resonance occurs at the radius where — (uj — fhVt) 2 + tz 2 = holds. After obtaining 
C£±n by solving this equation, we can obtain C^±h an d £ 2 n ±,n from equations (28) and (29), 
respectively. 

It is noted that in the coupling through n — n — 1 in warped disks (n = 1), which is the 
case to be considered below, we have A Z) ± t fi = (see the next subsection) and £, n ± n = 0. 

2. 3. Expressions for Coupling Terms 

Detailed expressions for the coupling terms, i.e., the right-hand side of equations (27) 
- (29), are very complicated (see Appendix for details). Their expressions, however, can be 
simplified by using the approximations of H <C A <C r. In the case of n = n — 1, after a lengthy 
manipulation we have 



nQ 2 ± dr 2 

dA z ^±^ 



(30) 
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,±,71-1 
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ip,iz,n— 1 



dQ 2 \ x ^ 2 TT d ( » d\ . , . - x d 



r \ dr 



'dr 




A z + n -i = small. 



(31) 

(32) 
(33) 



Here, A Zj ± tTl -i is small in the sense that H(dA Zj ± jTlf -i)/dr is negligible compared with A-^n-i 
in magnitude. 

Similarly, in the case of n — n + 1, we have 
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r,±,n+l 



yi, d 



dr dr r 



dr 



tw 



Ato ± n+i = smaller than A r ± n+ i in magnitude 



(34) 
(35) 



A z ,±,n 



+1 



r Sr.n "I - 'isr.n 7 ~F So),n 

otr ar r 



W* 



(36) 



It is noted that in the above expressions, ^ is adopted for A +iJ ^, while ^* is adopted 
in the case of A_ ^, where the asterisk denotes the complex conjugate. 

2.4- Growth Rate 

By the feedback process through the intermediate oscillations, the original oscillation is 
amplified or damped, so the frequency u can no longer be real. Kato (2008b) showed that the 
growth rate —0Ji, where u\ is the imaginary part of u, can be written as 
W± 



2E 



(37) 



where E is the wave energy of the original oscillation, W± is the rate at which work is done 
on the original oscillation by the nonlinear resonant process, and ± denotes the cases of the 

^ int ^ int 

coupling through £, and , respectively. Here, W± are written as follows, 



W 4 



cc'.r*; 



^3 J p m (r)[(27t) 3 ' 2 mrH] 



X 



pint A* _i_ pint A * 



dV, 



_pint 



n z,+,n 



dr, 



(38) 



W- 



-hit jW\ ^.-H' -int 



^s|p o(r)[(27r) 3 / 2 n!ri7] 



x 



pint J* 1 pint 4* _ pint 4* 

,n r,— ,h ' ^tp,—,h ip,—,h ' ~^z,—,h z, 



dr, 



(39) 



where Uq is the frequency of the original oscillation before the mode couplings, and the wave 
energy, E, is given by 



E = -u I Po C[oo-t(u-V)}^dV 
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(27T) 



3/2 



2 

where -E n is a dimensionless quantity given by 



u 2 (r*Hp 00 ) c E n , (40) 



J {rHp 00 ) c u V r c r c J 

with subscript c denoting the values at the resonant radius, where — (uq — rhVt) 2 + k 2 = 0. 

3. Results 

3.1. Original Oscillations 

As the original oscillations, we consider the axisymmetric g-mode oscillations trapped 
around the radius of the maximum epicyclic frequency, ft max . This type of trapped g-mode 
oscillations has been examined in previous studies (Okazaki et al. 1987; see also Kato 2001, 
Ferreira & Ogilvie 2008). For this purpose we solve equation (11) with m = and n — 1. In 
what follows, we assume that the sound speed (c s , in units of c) is constant. 

As mentioned previously, we adopt the radius of ISCO as the inner boundary radius. As 
the outer boundary radius, we take a radius in the evanescence region of the original oscillations. 
The amplitude of the oscillations is assumed to vanish at both inner and outer boundaries. 

Figure 2 shows the original oscillations for c s = 10 -3 . The upper panels shows the prop- 
agation diagram for the lowest two eigenmodes, the fundamental mode and the first overtone, 
although the two horizontal lines showing eigenfrequencies are indistinguishable on the figure. 
The left panel is for spin parameter a = 0.5, while the right panel is for a = 0.9. Note that the 
eigenfrequency is higher for a more rapidly rotating system, because the epicyclic frequency, k, 
increases with increasing spin parameter a. The eigenfrequency is also higher for a lower sound 
speed. The lower panels show the eigenf unctions of the fundamental modes (left panel) and 
the first overtones (right panel). The eigenmodes for the two different values of a look similar, 
but as a increases, the propagation region is shifted closer to the central object. 

3.2. Intermediate Oscillations 

To obtain the intermediate oscillations, we solve equation (30) with the nonlinear cou- 
pling terms. As discussed in the previous section, the azimuthal mode number of the interme- 
diate oscillation is given by rh = m ± 1 = ±1, and the vertical mode number n = ra±l = 0or 
2. By analysing the resonance condition for modes with Re(u) = ujq ~ ft max) we find that for 
in — — 1 and h = or 2, there is no resonance radius inside the propagation region of the original 
oscillation, so that there is no mode excitation in these cases. Therefore, for the intermediate 
oscillations, we have two cases of (rh,n) = (1,0) and (1,2) to study. 

In the case of (m,n) = (1,0), the intermediate oscillations are p-mode oscillations, and 
their propagation regions are those given by Uq < Q — k and > Q + k. The intermediate 
oscillations propagating in the region of ujq < Q — n have a resonance at the outer edge of the 
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Fig. 2. Eigenmodes and propagation diagrams of the original oscillations for c s = 10~ 3 . The up- 
per panels arc propagation diagrams showing the lowest two eigenmodes for a = 0.5 (left panel) 
and a = 0.9 (right panel). In the lower panels, the eigenfunctions are shown for the fun- 
damental modes (left panel) and the first overtones (right panel) of the original oscillations. 

region, and the resonant radius is inside or close to the propagation region of the trapped original 
oscillations, since u ~ ft max - Hence, we consider the intermediate oscillations propagating in the 
region of ujq < Q — k. The boundary conditions to be imposed to the intermediate oscillations are 
thus vanishing Lagrangian perturbation of pressure at the radius of ISCO and zero amplitude 
at the same radius where the outer boundary condition is imposed to the original oscillations 
(see subsequent paragraphs for implications of the latter boundary condition). 

In the case of (m,n) = (1,2), the interesting intermediate oscillations are g-mode oscilla- 
tions and their propagation region is Q — k < u>q < Q + k. At the inner edge of the propagation 
region, we have a resonance of ujq = Q — k and the resonant radius is inside or close to the 
trapped region of the original oscillations. When we impose the outer-boundary condition for 
the intermediate oscillations, caution is necessary, since the corotation radius exists inside the 
propagation region. 

At the corotation resonance, g-mode oscillations are damped (Kato 2003, Li et al. 
2003). Hence, this effect should be taken into account. However, a detailed consideration of 
the damping effect at the corotation resonance is complicated. Thus, a practical way is to 
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Fig. 3. Eigcnfunctions of the fundamental modes of original oscillations and the intermediate os- 
cillations obtained for Wq = 0.01. Panel (a) is for c s = 10 -2 and a = 0.5, panel (b) is for 
c s = 10 -3 and a = 0.5, and panel (c) for c s = 10~ 3 and a = 0.9. The black curve is the eigen- 
function of the original oscillation, the red and blue curves are respectively the real and imaginary 
parts of the intermediate oscillation. The vertical line denotes the resonance radius of ujq = O — k. 

introduce a damping of intermediate oscillations at the rate of (3Q with (3 being the damping 
parameter, as Ferreira & Ogilvie (2008) did. 

It is noted that unlike the case of (m,n) = (1,2), the effect of corotation resonance on the 
wave behavior will be unimportant in the case of (m,n) = (1,0) for the following reasons. In the 
case of p-mode oscillations, the corotation resonance can amplify the oscillations, different from 
the case of g-mode oscillations (Tsang and Lai 2009). Since the corotation radius is inside the 
evanescent region of the p-mode oscillations, the growth rate by corotation resonance depends on 
the efficiency of penetration of the oscillations through the evanescent region. In geometrically 
thin disks, the efficiency of tunnelling through the evanescent region is low (remember why 
the Papaloizou-Pringle instability does not have large growth rate in geometrically thin disks). 
Furthermore, in the case of m — 1, the evanescent region is wide compared with other cases of 
fh>2 (e.g., see Fig.4 of Tsang and Lai (2009)). Indeed, Tsang and Lai (2009) showed that the 
growth rate in the case of one-armed oscillations (m = 1) is much smaller than that in other 
cases. Thus, the effect of corotation resonance in the case of (m,n) = (1,0) can be neglected in 
our present problem. 

Considering these situations, we focus our attention in this paper only on the case of 
(m,n) = (1,0), with no consideration of the effects of corotation resonance, since in this case we 
can study the resonant excitation processes in an idealized form without introducing ambiguities 
related to the corotation resonance. 

Figure 3 shows the eigenfunctions of intermediate oscillations with (m, n) = (1,0) for 
c s = 10~ 3 and Wq = 0.01. Three cases of different sets of (c s , a) are shown. From the left panel, 
the adopted set is (c s ,a) = (10~ 2 , 0.5), (10 -3 , 0.5), and (10~ 3 , 0.9). In each panel, the real 



10 



5*10 



5x10 



0.01 



*- 5x10 J 0.01 0.02 



0.05 0.1 



0.1 



0.2 



0.5 



Fig. 4. Variation of the growth rates of the oscillations with the disk structure and spin parameter. The 
left panel shows the dependence of the growth rate on the sound speed for a fixed warp amplitude at 
the inner radius, Wq = 0.01, for a = 0.5 and a = 0.9, and the middle panel shows the dependence on 
Wo for a fixed sound speed c s = 10~ 3 , also for a = 0.5 and a = 0.9, while the right panel shows the 
dependence on the spin parameter, a, for a fixed sound speed c s = 10^ 3 , with Wo = 0.01 and Wq = 0.1 

and imaginary parts of the intermediate oscillation are denoted by the red and blue curves, 
respectively. The corresponding original oscillation (the fundamental mode) is also shown by 
the black curve in each panel. The vertical line on the y-axis denotes the resonance radius. 
Note that as c s decreases or a increases, the intermediate oscillations oscillate with shorter 
radial wavelength and confined to a narrower region around the resonance radius). Also note 
that the amplitude of the intermediate oscillations depend little on c s and a (see discussions in 
the subsequent subsection). 

3.3. Growth Rates 



To calculate the growth rate, —oj- u we need iterative processes. We first assume a trial 
value of —u)i, and solve the wave equation (30) to obtain £™+,ri an d Then, by using 

equation (37), we calculate the growth rate, If the resulting value of — U\ is not equal to 
the assumed value, we assume another trial value of — uj\ and repeat this procedure until both 
values coincide. Note that in the present case of (m,n) = (0,1) and (rh,n) = (1,0), the work 
integral, W + , given by equation (38) is written explicitly in the form of 1 
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ar \ uo — U 



uj — \l r 



2Vt 



d 2 W 
dr 2 



dr. 



(42) 



where and hereafter £ r , and £ z are used to abbreviate 9ft(£™+ ), £ r ,i, and £ 2j i, respectively 
To examine the parameter dependence of the growth rate, however, the following ex- 
pression for W+ is more perspective than equation (42). Let us multiply the complex conju- 



1 We sould not confuse W+ with the warp amplitude W. 
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gates of equations (27) and (28) by ^ nt and respectively, and sum them. Then, we have 
(■^A* + expressed in terms of £ (and £ ) alone. Taking the imaginary part of the 

expression, we have a simple expression for W + , which is 



W + = Y I Poo(r)[(27r)^rH] 



x2ui 



dr. (43) 



I, . P>\ l\ tint \i i |tmt|z\ ;o I tint tint* tint* tint \ 
^0 — sfHIsr I + 1^ I J - " Z Ur V ~~ *r ^ j 

If this expression for W+ is used, equation (37), which can be used to determine the growth 
rate is written as 



f J ' Pm {r)[{2^rH] 



(n - ^XICT + ICT) + - c% nt ; 



dr = E. (44) 



In this equation tu; does not appear explicitly. However, Ui is implicitly involved in £ , since 
£ depends on u\ through equations (27) and (28). Equation (44) is further simplified if A v 
can be neglected compared with A r , which is really allowed if the radial variations of £ r and ^ 
are rapid [see equations (31) and (32)]. In this case, we can adopt — (loq — ^)^™ t +22f2^ nt ~ 0, 
and equation (44) is reduced to 

^ |p o(r)[(27r) 3 / 2 rF] x (n - Wb)|gf*| a dr = E. (45) 

Calculations of the growth rate have been performed by using equation (42) for some pa- 
rameters characterizing the disk structure (warp amplitude at the inner radius Wq and acoustic 
speed c s ) and spin parameter, a, of the central source. Results show, as mentioned before, that 
the fundamental mode of the trapped oscillations is excited. Figure 4 shows the parameter 
dependence of the growth rate of the fundamental modes. The left panel shows the dependence 
on sound speed c s for spin parameter a = 0.5 and a = 0.9 with a fixed Wq = 0.01. This c s - 
dependence of the growth rate is mathematically understandable, if we consider equation (44) 
or (45). A smaller value of c s gives rise to shorter radial wavelengths of the original oscilla- 
tions [see equation (11)] as well as of the warp amplitude oscillations [see equation (16)]. This 
means that the coupling terms between the original oscillation and the warp, A, rapidly vary 
in the radial direction with large amplitude when c s is small, since the coupling has terms of 
the derivative of £ and £ with respect to r. Such large amplitude variations of the coupling 
terms give rise to large amplitude variations of the intermediate oscillations [see equation (30)], 
if the coefficients of ^ nt on the left-hand side of equation (30) are fixed. In other words, if 
the coefficients are fixed, condition (44) [or (45)] cannot be satisfied in the case of small c s . 
For condition (44) [or (45)] to be satisfied in the case of a smaller c s , the coefficients of ^ nt 
in equation (30) must increase in the case of a smaller c s so that the resulting ^ nt decreases 
and condition (44) [or (45)] can be satisfied. This is realized by increasing the imaginary part 
of LO{. Notice, for example, that the imaginary part of the left-hand side of equation (30) 
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involving £™r * s ~~ ^1(^0 — ^)£™r • I n summary, in the case of a smaller c s , the value of 
must increase so that the amplitude of the intermediate oscillations is roughly independent of 
c s and the condition (44) [or (45)] is satisfied. Numerical results shown in panels (a) and (b) 
of figure 3 are consistent with this argument. That is, comparisons of the results in panels 
(a) and (b) show that the radial wavelength of the intermediate oscillations becomes shorter 
as c s decreases. The amplitudes of the intermediate oscillations, however, are roughly on the 
same order in panels (a) and (b), which is realized by an increase of \ui\ in panel (b) compared 
with in panel (a). The result that the growth rate increases with a decrease of c s qualitatively 
agrees with that of Ferreira & Ogilvie (2008), but quantitatively our result slightly differs from 
Ferreira & Ogilvie's. For instance, for c s = 10 -3 , Wq = 0.01, and a = 0.5, our model provides 
I u>i |= 3.3 x 10~ 4 , while | Ui \= 4.2 x 10~ 4 in Ferreira & Ogilvie (2008). This difference is likely 
to result from the difference between our analytical approach and their numerical one. 

The middle panel of figure 4 shows that the growth rate increases with increasing warp 
amplitude. This is natural, since the oscillations are excited by having the interaction with 
the warped disk. Ferreira & Ogilvie (2008) discussed that the growth rate should grow with 
(dW/dr) 2 oc Wq for small Wo, because the excitation consists of two processes: i) the interaction 
giving rise to the intermediate oscillation and ii) feeding back to the original oscillation. In fact, 
their numerical results, which are for Wo < 0.02, agree with this argument. Our model, however, 
provides significantly flatter dependance of the growth rate on the warp amplitude than that 
obtained by Ferreira & Ogilvie (2008). This may be explained as follows: From the real and 
imaginary parts of equation (30), we see that if Ui was fixed, both the real and imaginary parts 
of C n +,n> i- e -' irr an d CT) w °uld be proportional to the coupling terms, which are proportional 
to the warp amplitude, W . Then the work integral, W+, (and the growth rate) would be 
proportional to Wq. As uj\ increases with increasing Wq, however, the terms proportional to uj\ 
in the coefficients of the inhomogeneous wave equations for £™* and [i.e., the terms coming 
from [{ujo — fhVt) 2 + K 2 ]£* nt in wave equation (30)] would become non-negligible in determining 
and This makes W+ (and the growth rate) deviate from W+ oc Wq. Instead, it leads 
to a dependence flatter than Wq as shown in the middle panel of figure 4. 

From the right panel of figure 4, we can see the dependence of the growth rate on the 
spin parameter of the black hole for a fixed sound speed (c s = 10~ 3 ). The growth rate increases 
with increasing spin parameter. This spin dependence of the growth rate will be due to a 
similar reason as that in the case of c s -dependence. That is, as the spin parameter increases, 
the warp rapidly oscillates in the radial direction in the inner region of the disk, since the 
difference between Q and Qj_ is large there [see equation (16)]. Let us assume tentatively that 
U{ is unchanged from the value of the case of a small spin parameter. Then, the intermediate 
oscillations have a large amplitude with a short wavelength, since the coupling terms have 
radial derivative of £^ [see equation (31) and (32)]. This means that the condition (44) [or 
(45)] cannot be satisfied. For the condition to be satisfied in the case of a large spin parameter, 
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fct^i ] must increase so that the amplitude of the intermediate oscillations decreases until the 

^ hit 

magnitude of £ becomes comparable with that of the case of a smaller spin parameter value. 
This will be a reason why the growth rate increases with an increase of the spin. The result 
regarding the a-dependence of growth rate qualitatively agrees with those obtained by Ferreira 
& Ogilvie (2008) for the fundamental g-mode oscillations. Note that for a = case, we found 
that there is no excitation of the fundamental g-mode oscillation. 

4. Summary and Discussion 

We have numerically studied the excitation of trapped g-mode oscillations in warped 
disks around black holes, based on a scenario proposed by Kato (2004; 2008a; 2008b). We 
first obtained trapped g-mode oscillations with eigenfrequencies close to the maximum of the 
horizontal epicyclic frequency. Then, we examined whether these modes are excited via the 
resonant coupling with the warp of the disk. We have found that the fundamental modes of the 
trapped g-modes are excited except in the case of non-rotating black holes. We have performed 
calculations of the growth rate of the fundamental mode for some parameter values such as 
the warp amplitude in the inner radius, the isothermal sound speed in the disk and the black 
hole spin parameter. We have found that the growth rate increases as the warp amplitude and 
spin parameter increases or the sound speed decreases. The dependence of the growth rate on 
these parameters qualitatively agrees with that obtained by Ferreira & Ogilvie (2008) for the 
fundamental mode, but quantitatively there are significant differences between the two. We 
suspect that these differences are due to differences of the model considered. 

Our results presented in this paper show that the excitation of trapped g-mode oscilla- 
tions by global disk deformation is rather sensitive to the disk structure, because the resonant 
radius given by cu = f2 — k and a characteristic radius of the oscillation where £ r . jn or d£ r>n /dr 
has a peak are close. We also suppose that the damping (or amplification) of the intermediate 
oscillations at the corotation resonance is not an essential ingredient for the excitation of disk 
oscillations in deformed disks. In the next paragraphs, we discuss the excitation mechanism 
from a viewpoint of energy flow among oscillations. 

Energy flow between original and intermediate oscillations occurs at two aspects of the 
resonance. The first occurs at the stage where the oscillation resulting from nonlinear coupling 
between the original oscillation and the warp acts as a forcing source on the intermediate 
oscillation. The second occurs at the feedback process from the intermediate oscillation to the 
original oscillation where the intermediate oscillation acts as a source to the original oscillation 
after coupling with the warp. The net energy input rate to the original oscillation through 
these coupling processes is given by equation (37). The results derived from this equation show 
that in a simplified situation the signs of wave energy of original and intermediate oscillations 
must be opposite for excitation of the original oscillation (e.g., Kato 2008a). This suggests that 
a physical cause of resonant excitation of disk oscillations in deformed disks is an interaction 
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of two waves with different signs of the wave energy at a resonant radius (Kato 2008a 2008b: 
Ferreira & Ogilvie 2008). In the present case, the original oscillation is axisymmetric and has 
positive energy, i.e., E > 0. Hence, a necessary condition for wave amplification to occur is a 
negative wave energy of the intermediate oscillation. The major process of excitation is thus 
a positive energy flow from the intermediate oscillation with negative energy to the original 
oscillation with positive energy at the resonant region. By this energy flow, both original and 
intermediate oscillations grow. 

The main role of the warp will be a kind of catalyzer for energy transport. In equa- 
tion (37), energy exchange between the oscillation and the disk in the resonant region will be 
involved, but it will be subsidiary in the present excitation mechanism. If the intermediate 
oscillation has an interaction with environment and its negative wave energy is transferred to 
the environment, the oscillation has a tendency to be damped. However, a negative energy will 
be transported to the intermediate oscillation from the original one to replenish the loss. Then, 
the original oscillation will be excited more strongly than that in the case of absence of such 
negative energy loss from the intermediate oscillation. This is the case considered by Ferreira 
& Ogilvie (2008). 

Here, a brief comment is necessary on the procedure adopted in this paper to evaluate 
the growth rate. Our procedure is based on the fact that the growth rate can be calculated 
by using W± given by equations (38) and (39). With these equations, we need not to solve an 
inhomogeneous wave equation for the original oscillations. In other words, the inhomogeneous 
wave equation to be solved by taking into account the resonant effects is only that of inter- 
mediate oscillations, as we have done in this paper. This is justified by comparing terms in 
equations (38) [or equation (39)] as follows. 

Let us, for example, consider 3(^ nt v4*) as a typical term in the integrand of equation 
(38) (the subscripts + and n are omitted for simplicity). The term consists of ^A^ + ^A^. 
If the effect of resonant processes on the original oscillation is neglected, A* is purely imaginary, 
i.e., A* r = [see equations (31) or (34)]. Then, 3(|j. nt A;) = ^A* r -. This is nothing but the 
case we have treated in this paper. If the effect of resonant processes on the original oscillation 
is taken into account, however, a real part of A* appears, i.e., ^A* r ^ 0. Then, the problem 
to be considered here is whether the term ^'I* r is of importance in evaluating W + . The term 
is, however, found to be smaller than C^A* ; in magnitude by the following reasons. 

A* T is non-zero, but smaller than A* ri in magnitude, since it is a small correction term of 
A* introduced by the resonant processes. Hence, the term £™fJl* r is negligible compared with 
^A*- v unless £™* is much larger than £™ in magnitude. The results in this paper show that 
£™* and ^ are of the same order in the resonant region (see figures 3). Therefore, ^A* r can 
be safely neglected compared with ^A*^ In conclusion, when we calculate the growth rate 
by using W±, the inhomogeneous wave equation that should be considered is only that for the 
intermediate oscillations, as we have done in this paper. 
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The results of this paper as well as those of Ferreira & Ogilvie (2008) confirm that 
the resonant excitation process of disk oscillations in deformed disks work to excite trapped 
oscillations. The excitation process is thus one of the most promising causes of the observed 
high frequency QPOs. The trapped g-mode oscillations themselves, however, cannot account 
for the observed characteristic that high frequency QPOs often appear in pairs with frequency 
ratio close to 3:2. Furthermore, in the case of kHz QPOs in neutron-star low-mass X-ray 
sources, the frequencies of QPOs vary with time. To describe such observational characteristics 
by the present disk oscillation models, it might be necessary to relax the restriction of trapping, 
i.e., consideration of non-trapped oscillations might be necessary, as Kato (2004, 2008a) did. 
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Appendix 1. Calculation of Nonlinear Coupling Terms 

The coupling terms used for the calculation is written as 

PoC(£,£i) = PoA i ' + p A p , (Al) 

where 

P „A*=-i^^(V«, (A2) 
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Their detailed expressions for r-, tp-, and z-components are given in Appendix 1 of Kato (2004). 
For the original oscillation we have 

£ r = exp[i(ut - m(p)YHn£ r ,n, (A4) 

£ v = exp[i(ujt - mip^Hni^n, (A5) 

£ z = exp[i(ut - rrup)]U n -i^ n , (A6) 

and for warps we have 

C = exp(-^)^ 1 | r M/ ) (A7) 

£jf = expH^il^, (A8) 

e = ex P (-i^)%lf , (A9) 

where ^ = -iW{r)H, = -W(r)H, and £f = iW(r)r. 

A. 1.1. Expressions for A in Cylindrical Coordinates 

Their derivations from equation (A2) are relatively easy. Using 0(£ Ztn ) — n(H/\)0(f- r)n ) 
and the approximations of H <C A <C r, we have, after some manipulation, 

( £ h (A10) 
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(A14) 
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(A15) 



A. 1.2. Expressions for A in Cylindrical Coordinates 

Derivation of detailed expressions for A^, and are somewhat troublesome. Below 
we consider them separately. 
A. 1.2.1. Expression for A? 

From equation (A3) we have 



PoA p r 
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T-rr ~t~ Turnip i 
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(A16) 



<9r rdip dz r 

where in the case of the coupling through , T rr , T^r, T zr , and T w in equation (A16) are 
respectively given by 
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For coupling through ^™ t , every £ w * in equations (A17)-(A20) should be replaced with £ w . 

Then, taking only the leading order terms under 0(j- z>n ) = n(H/X)0(Zr,n) and the ap- 
proximation of H <C A <C r, we have 



A-,±,n-l 



4 P 



V dr dr dr dr 



f2 2 ff^ r,n ^ 
± dr dr 



r dr 



A. 1.2.2. Expression for 

From equation (A3), we have 
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(A21) 
(A22) 
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where in the case of the coupling through T r(p and T zip are given, respectively, by 
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dz rdip dz rdip 

As in the expressions for A?, every £ w * in equations (A24) and (A25) should be replaced with 
£ w in the case of coupling through £™\ 

Substituting equations (A24) and (A25) into equation (A23), and taking only the leading 
order terms as we did before, we have 
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(A26) 
(A27) 
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A. 1.2.3. Expression for A p z 

The z-component of equation (A3) is written as 



PoA p z 



dT rz or 



dT z 



dr rd(p dz 



(A28) 



where in the case of the coupling through £™\ T rz , T vz , and T zz are given respectively by 
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For coupling through £™\ every £ w * in equations (A29)-(A31) should be replaced with £ w . 
Then, with the same procedures taking only the leading order terms as we did for A? and A^, 
the following expressions are obtained: 
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Using equations (A10)-(A15) and equations (A21), (A22), (A26), (A27), (A32), and 
(A33) with 0(£ z>n ) = n(H / XjO^n) and H <C A <C r, we finally have equations (31) - (36) given 
in the text. 
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